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ABSTRACT 

In this paper, we present an optimal open-loop slew trajectory algorithm developed at GSFC for the so-called 
“Yardstick design” [1] of the James Webb Space Telescope (JWST). JWST is an orbiting infrared observatory featuring 
a lightweight, segmented primary mirror approximately 6 meters in diameter and a sunshield approximately the size of a 
tennis court. This large, flexible structure will have significant number of lightly damped, dominant flexible modes. 
With very stringent requirements on pointing accuracy and image quality, it is important that slewing be done within the 
required time constraint and with minimal induced vibration in order to maximize observing efficiency. With reaction 
wheels as control actuators, initial wheel speeds as well as individual wheel torque and momentum limits become 
dominant constraints in slew performance. These constraints must be taken into account when performing slews to 
ensure that unexpected reaction wheel saturation does not occur, since such saturation leads to control failure in 
accurately tracking commanded motion and produces high frequency torque components capable of exciting structural 
modes. A minimum-time constraint is also included and coupled with reaction wheel limit constraints in the 
optimization to minimize both the effect of the control torque on the flexible body motion and the maneuver time. The 
optimization is on slew command parameters, such as maximum slew velocity and acceleration, for a given redundant 
reaction wheel configuration and is based on the dynamic interaction between the spacecraft and reaction wheel motion. 
Analytical development of the slew algorithm to generate desired slew position, rate, and acceleration profiles to 
command a feedback/feed forward control system is described. High-fidelity simulation and experimental results are 
presented to show that the developed slew law achieves the objectives. 
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INTRODUCTION 

JWST is an orbiting infrared observatory scheduled to be launched in 2011. Its deployed structure consists of a large 
flexible sunshield, on one side of which is a telescope with a lightweight 6.5-meter deployed primary mirror of 18 
segments and a suite of four science instruments, and on the other side is a spacecraft bus where the reaction wheels are 
mounted. For more details of the JWST configuration see [1]. 

One of the key performance metrics is the image quality that is driven by wavefront error (WFE) and image motion. 
The current allocation for WFE and image motion due to structural dynamics is 14 nanometers (RMS) and 4 milli-arc 
seconds (RMS), respectively. These stringent requirements dictate that dynamic disturbances induced internally be 
minimized to prevent any significant vibratory motion due to lightly damped system flexible modes that can exceed the 
allocation for required image quality. For JWST, the primary source of internal disturbance is from the reaction wheel 
assemblies. In addition to static and dynamic imbalances that continuously produce unwanted forces and torques 
proportional to the square of the wheel spin speed, [2,3], the reaction wheel torque and momentum limits when 
encountered with large torque commands will result in discontinuity in the produced torque profile and its derivative, 
and in turn, create significant vibratory motion in the system. This situation potentially can happen during slew when 


reaction wheels undergo a large swing in momentum to affect the motion of the spacecraft whose mass is much greater 
than that of the reaction wheel. Avoiding reaction wheel limits during slews not only maintains the internal disturbances 
within the image performance requirements, but also allows for rapid retargeting and better tracking capabilities by 
preventing control failure in tracking commanded motion during slews. 

This paper describes a slew trajectory algorithm developed to minimize both the effect of the control torque produced 
by reaction wheels on the flexible body motion and the maneuver time. With reaction wheels providing control 
actuation, this algorithm produces a desired slew trajectory in terms of position, rate, and acceleration profiles as 
functions of reaction wheel parameters, which are initial wheel speeds, and wheel torque and momentum limits of all 
the reaction wheels in the system. These reaction wheel parameters are taken into account in the algorithm to ensure that 
reaction wheel limits can only be tangentially reached during the slew so that no discontinuity in actual torque will 
occur. These kinematics profiles are used to provide commands to any stable feedback/feedforward control system in 
order to perform the maneuver. Note that this paper describes the slew trajectory algorithm only as an open-loop system 
and does not consider it in a closed-loop control design. 


DERIVATION 

There are two steps in the development of the algorithm. The first step involves selecting the basic slew trajectory, 
which defines and parameterizes the functional form of the desired slew position, rate, and acceleration profiles. This 
step, while using an approach similar to that described in [4], extends it by taking into account the reaction wheel limits 
in obtaining the desired motion profiles. The second step optimizes the slew command parameters, such as maximum 
slew rate and acceleration defined in the first step, for a given redundant reaction wheel configuration with torque and 
momentum limits, based on the dynamic interaction between the spacecraft and reaction wheel motion. 

In [4], a near optimal maneuver trajectory algorithm is given for the case with no reaction wheel constraints. This 
algorithm specifies, in a parametric form, a smooth function that basically approximates the sign function of a bang- 
bang law, which is a well-known form for minimum-time controllers. The trajectory motion in this case is continuously 
maintained at maximum acceleration over a long portion of the time during the slew period. The duration of this 
acceleration time portion is selected in a trade-off between minimum maneuver time and the degree of smoothness in 
the trajectory. 

In the presence of reaction wheel momentum and torque limits, arbitrarily long and continuously accelerating motions 
cannot be supported, implying that bang-bang controllers and the approach given in [4] are not applicable. With limits 
imposed on the reaction wheel momentum, the best that can be achieved in terms of minimizing maneuver time is to 
maintain spacecraft motion at its maximum rate allowable by these limits, over as long a portion of the slew time as 
possible. 

Maximizing this coasting time period means minimizing the time it takes to reach the maximum rate at the start of the 
slew, and the time it takes to go back to zero rate at the end of the slew. This requires that the maximum allowable 
acceleration and deceleration, which are bounded by the reaction wheel torque limits and initial wheel momentum, be 
used in getting to and from the maximum rate. 

Thus, the general acceleration profile appropriate for this case starts with an impulse whose amplitude is derived from 
reaction wheel torque limits, and whose direction is the direction of the desired motion. Following the impulse, the 
profile remains at zero over the entire slew duration, and ends with an impulse equal to the one applied at the start but in 
the opposite direction. The resulting rate profile is simply a rectangular curve, whose amplitude is derived from the 
reaction wheel momentum limits. 

In order to minimize the structural mode excitation, a smooth version of the acceleration impulsive profile described 
above is selected, and is chosen so that it and its first derivative are continuous and expressible in low order 
polynomials. The basic polynomials used in [4] are employed to provide smooth slew motion including at the two end 
points. 
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Derivation of Slew Profile 


To achieve a smoothly varying torque, the desired acceleration profile generated by the slew trajectory algorithm has 
the following functional form: 

E.l a(t) = a d e~f(t,At a ,At c ) 

at 

Where Od is the maximum acceleration amplitude, e is a constant vector specifying the slew rotation axis, and/( ) is a 
positive scalar time-function representing the desired spacecraft rate, parameterized by (At a Atc), and whose first order 
derivative is given as: 


E.2 


4/«. &.,&,)=< 

at 



3-2 


3-2 


v At « JJ 


V At a JJ 


3-2 


J t - t W 


V Al a J 


\ 2 ( r ts _ t w 


V At a JJ 


3-2 


t 0 <t <tj — tg + A t Q 

tj < t < t 2 = t, + A t a 
t 2 < t < t 3 — 1 2 A t c 
t 3 <t<t 4 =t 3 + A t a 

t 4 <t<t s =t 4 + A t a 

t>t 5 


The time periods At a and At c can be expressed in term of the maximum acceleration magnitude £& , the maximum rate 
magnitude tiJi , and the slew angle A# as: 
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and the total slew period is: 
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It follows that the desired rate profile has the following form: 

E.5 co(t) = a> d ef(t,At a ,At c ) 


The profile of the slew angle Zl#can be directly obtained by integrating E.5 to give: 

T 

A 0(T) = co i e \f(i,At a ,At c )di 
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Figure 1 shows the normalized slew acceleration, rate, and angle profiles described by El, E5, and E6, when both a d 
and co d are set to 1. Also, in these plots, the time periods are selected to be At a =1 second and At c , =2 seconds. 
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Figure 1 - Idealized Trajectory Profile 

Since the dead band time At a cannot be negative, the maximum rate magnitude Q) d must also satisfy the following 
condition: 

. la, A<9 
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which can be given in terms of At a by using the first relation in E.3 for C& , i.e.: 
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For later reference in the derivation of the maximum acceleration magnitude, given here are the relations between Qd, 
COd, At a , and the slew angle A6( At a ), at time At a after the start of the slew, when the desired acceleration profile 
reaches its positive maximum amplitude <3^, and the desired rate profile reaches half of its maximum value CQd : 
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Derivation of the maximum rate magnitude 

The basic equation describing the system momentum of the spacecraft and reaction wheels in the inertial frame is: 


E.9 


H 0 = R T (t)(Alh w (t) + I s m) 


Where h w is the wheel momentum vector in wheel frame, and A * represents the constant transformation from the wheel 

frame to the spacecraft frame. In a four-wheel configuration, h w is a 4-component vector, and A^ is a 3x4 matrix whose 
columns are reaction wheel rotational axes given in the spacecraft frame and can be obtained directly from the reaction 
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wheel configuration. The left inverse of this matrix, denoted by A £ , is a 4x3 matrix and typically given to minimize the 

wheel speed magnitude. Ho and I s are, respectively, the initial system momentum in the inertial frame, and the system 
inertia matrix in the spacecraft frame. Ho is assumed here to be a constant vector, i.e. there is no external torque, and 
consists only of wheel momentum. R(t) is the direction cosine matrix representing the transformation from the inertial 
frame to the spacecraft frame corresponding to the desired rate profile CtXf). 

Using the minimum mean squared solution, A£ , for the inverse of A w , the wheel momentum vector of minimum 

magnitude given in the reaction wheel frame, as a function of the initial system momentum and the desired motion 
profile can be written as: 

E.10 h w = AZ{R(t)H 0 -I s m) 


For a redundant wheel configuration, where there are more than 3 wheels in operation for attitude control, if the null 
vector in the wheel frame is non-zero initially, then it must be added to E.10 and must be carried in all the subsequent 
steps in computing the maximum rate magnitude. 

With Equation E.3, E.10 can be expressed in terms of slew parameters as: 
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Let ±h lim and ±T hm , h Vlm >0 and Zf lm >0, denote the momentum limits and torque limits, respectively, of each reaction 
wheel. These limits correspond to the four comers of a rectangular torque-speed curve. Although a typical torque-speed 
curve may have regions of decreasing torque as the wheel approaches saturation speeds, only the rectangular torque- 
SDeed curve is considered here. 


The momentum limits when applied to a reaction wheel leads to the following inequality: 

E - 12 ~Kn * K i =(AZ) i {R(t)H 0 -I s eCO i f(t)) < 

Where the index i indicates the i“ wheel, and ( )j denotes the i* row of the body-to- wheel transformation matrix. 

E.12 leads to the following upper and lower bounds of the spacecraft rate component along the 1 th wheel axis, due to the 
momentum limits of the I th wheel: 

£■13 -K. + S (AZ)lMf<f) £ K, + (A"l «(*)». 

Global bounds over the slew period can be obtained by taking the appropriate maximum and minimum values of the 
time-dependent terms over the period, i.e. 

E 14 -*,„+Max{(4') < (A'llto/m S + ),*«».} 

The expression inside the above Min and Max operators taking over the slew time period represents the component of 
the system momentum that the i* wheel must be able to carry during the slew. As the spacecraft rotates, there are for 
each reaction wheel, two attitudes at which the system momentum on the wheel are maximum and minimum, which 
occur when the total momentum vector in the spacecraft frame is closest and farthest from the wheel’s rotational axis. In 
addition the initial and final slew attitudes may represent extreme points in the case where there is no extreme point 
during the slew. For the four wheel case, there are 10 of these rotation angles, and the exact solutions of these angles 
can be obtained since R( •) can be expressed in terms of position without being explicitly dependent on time. 

The computation of these attitudes is simplified by considering a coordinate frame, referred to as the slew frame, in 
which the z-axis coincides with the slew axis, e , and the x axis is defined such that the initial system momentum H 0 is 
in the x-z plane. In this frame, motion of the vector H 0 describes a cone about the z-axis as the slew progresses, and for 
a given wheel rotational axis vector, the minimum and maximum functions of E.14 can be computed straightforwardly. 
They can be shown to correspond to having Ho at either the end points or lying in the plane spanned by the z-axis and 
the I th wheel’s rotational axis. 

Since /(•) is a positive function with the two end points at zero value, the upper and lower bounds on the body rate 
given in E.14 must have one bound positive and the other negative. Thus, the maximum rate amplitude, denoted by 
G ) di , based on the i* wheel’s momentum limit is simply the larger of the two bounds (i.e. the positive bound). That is: 
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Also, the condition for a valid command slew is that the upper and lower bounds defined in E.14 must have opposite 
signs. This assures that the component of the system momentum on a reaction wheel cannot exceed h Vim during the slew. 

The sign of the denominator in E. 15 plays a very important role in determining which of the two bounds is the 
maximum rate amplitude, as it indicates the direction in which that wheel must be rotated to achieve the desired 
spacecraft motion in the absence of the initial momentum. In general, it has the same sign as that of the wheel 
momentum limit that will be reached during the slew. In the case where this quantity is identically zero for the 1 th wheel. 



which means that there is no contribution from this reaction wheel to the slew and thus, it can be removed from the 
computation of the slew maximum rate and acceleration magnitudes. 


The maximum rate amplitude (O d that satisfies the momentum limits of all the reaction wheels in the system, is the 
smallest of all the CO d . : 

E.16 co d = Minjfl^} 

In the special case, where the initial system momentum H 0 is along the slew axis e , Equation E.13 for all the reaction 
wheels, reduces to: 
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Which has constant bounds so that the maximum rate magnitude is easily obtained in this case. 


Derivation of the maximum acceleration magnitude 

The wheel torque required to carry out the desired motion is obtained by taking the derivative of E.10, which gives: 
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Where the square brackets denote the matrix corresponding to the cross product operator, i.e. [a] b = a X b . 


With Equations E.l and E.3, E.17 can be rewritten in terms of slew parameters as: 
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Using the same convention as in the above section in expressing components related to a reaction wheel, the following 
inequality is obtained when applying the torque limits to the i* wheel: 

E.19 -r„ = «i. 

The index i appears in Ofa to indicate that the derived parameter in this case is the maximum allowable acceleration due 
to the 1 th wheel. Rearranging E.19 leads to the following upper and lower bounds on the component of the spacecraft 
acceleration along the i ^ wheel axis, due to the torque limits on that wheel: 
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At this point, an approach similar to what was used previously to derive the maximum rate magnitude can be used to 
obtain a legitimate maximum value for the acceleration profile. This involves obtaining the global bounds for the time- 
dependent term in E.20, which is the gyroscopic effect term, over the slew period, and taking the minimum of the norm 
of the two bounds as the maximum allowable acceleration by the i* wheel. The desired acceleration magnitude is 
simply the minimum of these maximum allowable acceleration values over all wheels, i.e. by simply setting f(t) and its 
time-derivative respectively to ¥2 and ±1, which are the values when the acceleration is at its maximum values, Equation 
E.20 can be manipulated to yield the following maximum wheel accelerations: 
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so that: 


a d = Min{a dl ] . 


This is an admissible solution because as the acceleration profile, which is a function of the derivative of j\t) given in 
Equation E.2, changes sign from positive to negative in the course of the slew, the norm of the smaller bound will 
ensure that the acceleration profile will not exceed torque limit of the i* wheel. By taking the minimum over these 
smaller bounds of all the reaction wheels, the acceleration profile will be within all the wheel torque limits. Although it 
produces a legitimate solution for the maximum acceleration magnitude, and is actually quite easy to implement, this 
approach in general results in sub-optimum acceleration profile. 


An approach, which yields an optimum solution for the maximum acceleration amplitude, derives the maximum 
acceleration values for the positive and negative regions of the acceleration profile separately, and selects the minimum 
of the two as the solution. This approach turns out to be more complicated to implement since it requires an iterative 
method for finding roots of an equation involving transcendental functions. 

The inequality given by Equation E.20 when applied to the positive region of the acceleration profile reduces to: 
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Since Equation E.21 is true for time between the starting time t 0 and t 0 + 2 At a , evaluating the above inequality at time 

At a when the acceleration profile reaches its positive maximum magnitude and the derivative o fj\t) equals to 1, leads 
to the following inequality: 
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This inequality contains two unknown parameters, 6fo(+) and At a , which, when obtained by solving the equality of 
E.22, are the optimum allowable acceleration and the shortest time to reach the maximum rate magnitude, satisfying the 
system torque constraints and relations in E.5 for the i A wheel in the positive acceleration region. Although it is possible 
for Equation E.22 to have more than one set of solutions, the set which is the desired solution has the largest £&(+) 
and smallest At a . 


With Equation E.8, the equality of E.22 can be rewritten in terms of just one variable, the slew angle evaluated at 
At a , for the I th wheel: 
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By using the slew coordinate frame defined in the previous section, E.23 can be simplified to the following form: 
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Where A0 a = is the slew angle at which the acceleration profile attains its positive maximum magnitude CCd and 

the rate profile attains half of its maximum magnitude Qfy and A , B , C, and <j> are constants corresponding to terms in 
Equation E.23, i.e. 
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Where S^*) and S 2 (0 are scalar functions given in terms of specified parameters, with Si assumed to have values 
between (-1,1), which limits the possible size of C. The terms of E.24 can be visualized as follows: B is the maximum 

possible acceleration that the i* wheel can provide along the slew rotation axis e . The term C sin(Aft + (f ) ) is the 

acceleration required to commutate the system momentum through the wheels as the slew progresses. The left-hand side 
of Equation E.24 is the remaining acceleration available to slew the body. 


Although Equation E.24 can not be solved analytically for A0 a , some simple iterative search methods can be devised for 
finding the smallest root, which corresponds to the desired (minimum slew time) solution. 


One adequate search technique to find the first root is to start at the top left of the region, which corresponds to the zero 
angle and acceleration (Z?+|Cj). The angle corresponding to this point using the left-hand side of E.24, i.e. 
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, should be within the bounded region and on the left of the first root. Starting with this angle, by 


stepping in an adequately small increment along the curve of the left-hand side of E.24, the interval, which contains the 
solution, can be determined. The following iteration method can be shown to yield the desired solution within this 
interval: 
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However, at any point in the above process, if an angle d i > 0 max , where 6 W is the angle defines the boundary of the 
first quadrant of the total slew angle, the maximum rate magnitude must be decreased, according to E.7, to give a 
smaller A that will permit the above process to continue. This rate reduction is simply to adjust the desired motion 
profile so that available reaction wheel torque in the system can allow the motion to accelerate to the maximum rate 
within the time constraint to have a dead band in the acceleration profile. 

The same method described above can be applied to the negative region of the acceleration profile to obtain the desired 
solution for the negative region. The equivalent equation to E.23 for the negative region of the acceleration profile is the 
following: 

WA") [c]fi(A0(A/,-AO)H„ 

Having obtained the maximum allowable acceleration values for each of the reaction wheels in the positive, <%/(+), 
and negative, Cfe(-), regions, the maximum acceleration magnitude is the minimum value taken over all these values, 
i.e. 
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In the special case, where the initial system momentum H 0 is along the slew axis e , the properly modified version of 
E.23 and E.27 for the i* wheel, which has non-zero term l s e , leads to the following solution: 



Where (Dd and (Xd are computed from the above formulation based on reaction wheel configuration, momentum and 
torque limits, and slew command, and completely define the desired slew trajectory. 


SIMULATION EXAMPLE 

A simple time domain simulation was created to demonstrate spacecraft slew performance using this trajectory 
algorithm. The simulation was implemented using rectangular integration with a 1 second time step. The system model 
included only rigid body dynamics of the spacecraft, with the following inertia matrix (from an early conceptual design 
for JWST) 
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The system model also included PD feedback control with a 1 Hz bandwidth, and feedforward control (to X h) to 
compensate for gyroscopic effects. The reaction wheel assembly in this model was a 4-wheel configuration. Each wheel 
was modeled using a rectangular torque-speed curve, with a 60 Nms momentum limit, a 0.3 Nm torque limit, and a 
rotor inertia of 0.0948 kg m 2 . The initial wheel momentum vector was (53.608, 41.695, 1 1.913, 5.956) Nms. The 
spacecraft was commanded to slew 359 degrees about the Y axis. The results are shown in figures 2 and 3. 

The total time to complete the slew was approximately 1800 seconds, as can be seen from the acceleration, velocity, and 
angle trajectories in figure 2. Comparing figure 2 to figure 1 shows that the salient features of the ideal trajectory are 
present, namely the smooth transitions between the segments of the trajectory, the smooth peak in the acceleration, and 
the flat coasting at maximum velocity. Figure 3 shows that hard saturation in both wheel torque and wheel momentum 

is avoided. The torque limit is approached in wheel 1 (at t = 500 seconds) and in wheel 2 (at t ^ 1500 seconds). The 

momentum limit is approached in wheel 2 (at t = 1150 seconds) and in wheel 4 (at t = 650 seconds). Thus, excitation of 
flexible modes would be mitigated. 
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Figures 2 - Slew Trajectories from Simulation 
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Figure 3 - Wheel Torque (left) and Wheel Momentum (right) Profiles from Simulation 













CONCLUSION 


An optimal open-loop slew trajectory algorithm is described. This algorithm enables minimum-time maneuvers of a 
spacecraft while avoiding the excitation of structural modes that would otherwise impact settling time. Furthermore, the 
algorithm does this without knowledge, or a model, of the flexible structural dynamics - all that is required is 
knowledge of the spacecraft rigid body mass properties and the reaction wheel torque and momentum limits. It has been 
implemented and tested via simulation to support the Yardstick design of the JWST. Additionally, experiments were 
performed with the MIT Origins Testbed, using a reduced-order form for a single-axis system, which confirmed the 
ability to track the commanded trajectory and to avoid excitation of unmodeled dynamics. 
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